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ABSTRACT 

A new Chandra observation of SNR 0506—68.0 (also called N23) reveals a complex, highly structured 
morphology in the low energy X-ray band and an isolated compact central object in the high energy 
band. Spectral analysis indicates that the X-ray emission overall is dominated by thermal gas whose 
composition is consistent with swept-up ambient material. There is a strong gradient in ambient density 
across the diameter of the remnant. Toward the southeast, near a prominent star cluster, the emitting 
density is 10-23 cm""^, while toward the northwest it has dropped to a value of only 1 cm~'^. The total 
extent of the X-ray remnant is 100" x 120" (24 pc x 29 pc for a distance of 50 kpc), somewhat larger 
than previously known. The remnant's age is estimated to be ~4600 yr. One part of the remnant shows 
evidence for enhanced O, Ne, and perhaps Mg abundances, which is interpreted as evidence for ejecta 
from a massive star core collapse supernova. The compact central object has a luminosity of a few times 
^q33 gj-gg g-i g^j^j obvious radio or optical counterpart. It does not show an extended nebula or pulsed 
emission as expected from a young energetic pulsar, but resembles the compact central objects seen in 
other core collapse SNe, such as Cas A. 

Subject headings: ISM: individual (LHA 120-N 23, SNR 0506-68.0) — shock waves — supernova 
remnants — X-rays: individual (CXOU J050552.3-680141) — X-rays: ISM 



an energetic pulsar or ejecta-dominated emission), while 
nearly all larger remnants are dominated by swept-up am- 
bient medium (Hughes et al. 1998). Of course, size alone is 
an imperfect predictor of evolutionary state, yet the pos- 
sibility that SNR 0506—68.0 is in a transitional stage of 
evolution (between ejecta-dominated and Sedov) makes it 
worth studying. The original goals of our observation were 
to find and categorize X-ray-emitting SN ejecta, search 
for a pulsar-powered synchrotron nebula or other com- 
pact X-ray star, and clarify the evolutionary state of SNR 
0506—68.0 through spatially resolved X-ray spectroscopy 
to better understand the nature of its progenitor star and 
its effects on the local environment. While this work was 
being completed, an independent report (Hayato et al. 
2006a, b) noting a spectrally hard. X-ray compact object 
in N23 came to our attention. In the following we present 
the results of our comprehensive exposition of the Chandra 
data in which we arrive at a different conclusion regarding 
the nature of the compact star. Throughout, we assume a 
distance of 50 kpc to the LMC. 

2. DATA REDUCTION AND ANALYSIS 

We observed SNR 0506-68.0 with Chandra on 2002 Dec 
29 using the ACIS-S detector in full-frame timed exposure 
mode (Obsid 2762). Standard cleaning procedures were 
applied to the data including filtering for times of high 
background that resulted in a useful live-time corrected 
exposure of 34744 s. Event pulse heights were corrected 
for time dependent gain effects and spectral response func- 
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1. INTRODUCTION 

The Chandra X-ray Observatory has been a revolution- 
ary tool for studying supernova remnants (SNRs) in the 
Large Magellanic Cloud (LMC). The telescope's superb 
spatial resolution has revealed newly processed metals in 
the interiors of most LMC SNRs and at least one new pul- 
sar wind nebula (see Weisskopf & Hughes 2005, and ref- 
erence therein). Studying the full sample of LMC SNRs 
is essential to gain a more complete understanding of the 
evolution of metal-rich ejecta from its synthesis in stars 
and SNe to its integration into the interstellar medium. 

Although relatively bright in the radio and X-ray bands 
(it is ranked ninth brightest in the soft X-ray band and 
eleventh brightest at 408 MHz; Mathewson et al. 1983), 
SNR 0506-68.0 (also called N23) has not been extensively 
studied. Chu & Kennicutt (1988) suggest that it had a 
massive star progenitor due to the high density of OB 
stars in the environment and its close proximity to H II 
regions and star clusters. There is a strong gradient in 
brightness across the remnant in the X-ray (Mathewson 
et al. 1983) and radio (Dickel & Milne 1998) bands. Its in- 
tegrated X-ray emission (Hughes et al. 1998) is dominated 
by soft thermal emission from swept-up ambient medium 
with the lowered abundances (^0.3 solar) typical of the 
LMC. 

According to Mathewson et al. (1983), SNR 0506-68.0 
is an intermediate sized remnant ('-~^11 pc in diameter). 
All smaller LMC remnants show evidence for youth (i.e.. 
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tions accounted for absorption due to contamination build- 
up on the ACIS filters. The optical field of the SNR is quite 
crowded and it was not possible to identify unambiguous 
counterparts to serendipitous Chandra sources and thereby 
improve the astrometric accuracy of the X-ray data. Nev- 
ertheless the observation is not affected by any known as- 
pect offset and therefore we expect the absolute positional 
uncertainty to be 0.6" (radius, at 90% confidence). 



0.5-2 keV 




2-6 keV 




05^ 06"^ 0^ 05"^ 45^ 

Fig. 1. — Soft (top) and hard (bottom) band Chandra X-ray 
images of SNR 0506—68.0. Spectral extraction regions are shown 
in green in the top panel. Contours delineate the faintest outer 
extent of the remnant from a smoothed image. 

2.1. Gaseous Remnant 

Figure 1 shows Chandra images of SNR 0506—68.0 in 
two broad energy bands (0.5-2 keV and 2-6 keV) along 
with contours indicating the faint outermost extent of the 
remnant from an adaptively smoothed image (not shown). 
Nearly all of the X-ray emission occurs below 2 keV ex- 
cept for a central compact object (see below) and some 
of the brightest portions of the limb. Those bright por- 
tions lie mostly toward the east and southeast and ap- 
pear highly structured and filamentary. The northwestern 
half of the remnant is much fainter and shows a hint of 



limb-brightening. The overall extent of the remnant is 
100" X 120", larger than previous estimates of the X-ray 
size (Mathewson et al. 1983; WiUiams et al. 1999). The 
integrated X-ray spectrum (Fig. 2) clearly indicates the 
characteristics of soft thermal emission and closely resem- 
bles the ASCA spectrum (Hughes et al. 1998). 

TABLE 1 

Spectral Model Fit Results for SNR Filaments 



Parameter 


Value and Uncertainty (1 a) 
East rim Soutlieast rim 


7Vh(LMC) (H atoms cm-^) 

kT (keV) 

logio (net/cm"^ s) 

Oxygen 

Neon 

Magnesium 
Iron 

nenHV/(47rD2) (cm"'^) 
xVd.o.f 


3.3l° J? X 10^1 
0.48 ±0.14 
10.6«- 
0.15 ±0.05 
0.23 ±0.09 
0.18 ±0.13 
0.13 ±0.06 
9.75 X 10'° 
94.8/56 


2.lt° * X 10^1 
0.77 ±0.4 
11.2 ± 0.3 
0.30 ±0.18 
0.32 ±0.09 
0.21 ±0.17 
0.20 ±0.15 
2.28 X 10'° 
42.4/49 


Two representative spectra from portions of the bright 



eastern rim (hereafter east rim and southeast rim, indi- 
cated by the two ellipses on Fig. 1) were extracted, re- 
grouped to contain at least 20 counts per channel, and an- 
alyzed with planar-shock nonequilibrium ionization (NEI) 
models (Hughes et al. 2000). The fits included Galactic ab- 
sorption with solar composition and column density fixed 
to the value 5.8 x 10^° atoms cm^^ from Dickey & Lock- 
man (1990) and LMC absorption with variable iVH(LMC) 
and metal abundances of 0.3 times solar. Numerical re- 
sults are given in Table 2.1 and the spectral data with 
best-fitted models are plotted in Fig. 2 (as blue and green 
curves). Only the abundances of O, Ne, Mg, and Fe in 
the emission model were allowed to vary freely, all other 
species had their abundances fixed to 0.3 times solar. The 
derived abundances clearly show that the bright rim emis- 
sion is dominated by gas with LMC composition. We also 
extracted the spectrum of the faint western limb as well 
(from the large rectangular region to the west of Fig. 1). 
All elemental abundances were fixed to 0.3 times solar, be- 
cause of the low statistics of this spectrum (~220 counts 
above background). Furthermore the uncertainties on 
the fitted thermodynamic quantities are large, so we only 
quote best fit values as indicators of the thermodynamic 
state: iVH(LMC) ~ 7 x lO^i H atoms cm'^, kT ~ 0.3 keV, 
logio(neVcm-3s) ~ 9.8, and nennV/iAnD^) - 4.5 x 10^° 
cm~^. The fit was statistically acceptable {x^ = 7.2 for 8 
degrees of freedom). 

Numerous spectra were extracted from elsewhere across 
the remnant including both filament and inter-filament re- 
gions, but in nearly all cases the gas composition was con- 
sistent with LMC abundances. Only one region showed 
any evidence for elevated abundances (the inverted "L"- 
shaped region in Fig. 1). This spectrum's relatively strong 
emission around 0.6 keV (see cyan curve in Fig. 2) clearly 
sets it apart from the rim spectra discussed above. Our 
spectral fit prefers high temperatures (>0.7 keV) and en- 
hanced abundances for O (>0.7) and Ne (>0.8), and per- 
haps Mg (>0.3). Note that these values represent the 
90% confidence lower bound on the parameters. Over this 
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range the Fe abundance is clearly constrained to be at 
or below the LMC value. For temperatures lower than 
0.7 keV, the abundance ratio of O to Fe (relative to so- 
lar) grows even larger and at 0.5 keV the ratio is >20. 
Although the precise values are somewhat indeterminant, 
we conclude with some confidence that the low Z elements 
are enhanced in this particular part of the remnant. 



2.2. Compact Star 

Close to the geometric center of the SNR an isolated 
X-ray point source (hereafter CXOU J050552. 3-680141) 
is clearly seen in the 2-6 keV band (Fig. 1) at a position of 
RA=05:05:52.32, Decl.=-68:01:41.2 (epoch J2000). The 
detection in this hard X-ray band is significant; we obtain 
31 net counts within a 2 pixel radius with an expectation 
value from background alone of 1.4 counts. The source 
is also detected in the soft band (0.5-2 keV) at approxi- 
mately the same significance level. The net (local back- 
ground subtracted) count rate of the source in the 0.5-6 
keV band is (2.8 ± 0.4) x lO'^ s'^. 




0.5 1 

E (keV) 



Fig. 2. — Spectra and best fit models for several positions in SNR 
0506—68.0. In order from the top the spectra are from the entire 
remnant (black), the east bright rim (blue), southeast rim (green), 
enhanced abundance region (cyan), and the central compact object 
(red). 



We simulated the Chandra point-spread-function (PSF) 
at the observed position of the compact object using Chan- 
dra Ray Tracer (ChaRT).*^ The observed structure and ra- 
dial profile of the star are entirely consistent with the PSF 
and we find no evidence for extended emission (for a close- 
up view of the point source see the insert box in Fig. 3). 

^ ChaRT is available on-line at: http://asc.harvard.edu/chart/index. 



TABLE 2 

Spectral Model Fit Results for Compact Source 





\ / Cl 1 1 1 O Tl H 

V ojiLlt; dliLl 


Parameter 


Uncertainty (1 o") 


JT UW til J-JcLVV iViUU-tli 




AfH(LMC) (H atoms cm~^) 




ap 




FE(lkeV) (photon s"^ cm''^ keV"^) 
xVd.o.f 




10.0/5 


Blackbody Model 




iVH(LMC) (H atoms cm^^) 


5.7t; j X 10^1 


kT (keV) 


U.aD_Q 20 


Norm (i?/lkm)V(-D/10kpc)2 


3.5l^:^ X 10-' 


xVd.o.f 


12.9/5 
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Fig. 3. — Broad band blue optical image from the MACHO 
project in the vicinity of the compact central X-ray source. Stars 
labeled A and B are closest to the X-ray source, whose error circle 
is shown in white. The X-ray data in the 2-6 keV band are shown 
in the small insert figure to the lower right. The slight E-W ex- 
tent of the X-ray source is entirely due to the Chandra PSF at this 
location. 

The spectrum of the compact source was extracted from 
a circular region 2 pixels (0.984") in radius. The main 
source of background comes from thermal emission emit- 
ted by the gaseous remnant itself, which accounts for half 
the extracted counts (~100 out of 193 total). Our spectral 
fit for the compact star, therefore, included a NEI thermal 
component with variable intensity, but fixed abundance 
(0.3 times solar), temperature and ionization timescale. 
The latter two quantities were determined by a spectral fit 
to data extracted from an annulus surrounding the point 
source covering 4-6 pixels in radius. Emission from the 
compact star was assumed to be either a blackbody or 
a power-law. Absorption, including Galactic and LMC 
components, was included as above. Table 2.2 shows the 
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results and the spectral data and model are plotted in 
Fig. 2 (in red). The red dashed curve shows the best fit 
power-law model for the compact star's emission. The lu- 
minosity of the compact source under the power law spec- 
tral model is Lx ~ 2 - 5 X lO^^ erg s^^ (0.5 2 kcV) and 
Ljf ~ 3 - 6 X 1033 erg s"! (2-6 keV). The blackbody fits 
yield similar values: ^ 0.8 - 2 x 10^3 erg s~^ (0.5-2 
keV) and ~ 3 - 6 X 1033 erg s"! (2-6 keV). 

Using FFT and epoch folding techniques we searched 
for a pulsed signal from the compact object. Nothing sig- 
nificant was found for frequencies of several mHz up to the 
Nyquist limit (~0.154 Hz), which was set by the intrinsic 
time resolution (^3.24 s) of our ACTS observation mode. 
We are unable to set interesting limits on the pulsed frac- 
tion because of the faintness of the source in comparison to 
the contaminating remnant flux even within a ~1" radius 
extraction region. 

Radio images (Dickel & Milne 1998) from the ATCA 
show no evidence for the compact source in any of the 
four wavelengths mapped (3.5 cm, 6 cm, 13 cm, and 20 
cm) to an estimated flux limit of ~2 mJy. These limits 
are a factor of ~20 below the radio flux of the pulsar wind 
nebula in SNR 0453-685 (Gacnsler ct al. 2003), whereas 
the X-ray source in SNR 0506—68.0 is only a factor of 
two fainter than the source in SNR 0453—685. Clearly 
CXOU J050552.3-680141 is relatively radio quiet. There 
is also no clear optical counterpart in the DSS, but the field 
is crowded and the flux limit is high. Wc obtained public 
domain images from the MACHO Project that allowed us 
to set a more sensitive limit. In Fig. 3 we show a portion 
of an image made by stacking 10 individual 300-s expo- 
sures with moderately good seeing. The combined image 
was registered to the sky using several hundred HST guide 
stars for a positional accuracy of <0.3" in each axis. A sub- 
set of these stars was also used to estimate a coarse photo- 
metric zero point for the image (accurate to roughly ±0.3 
mag). The nearest stars are 3.4" (star A) and 1.1" (star 
B) away from the X-ray compact object and are therefore 
nominally outside the Chandra error circle as indicated by 
the circle in the figure. The estimated magnitudes of these 
stars are ~17 (A) and ~19 (B). 

3. CONCLUSIONS AND SUMMARY 

There is a strong gradient in density from east to west 
across SNR 0506—68.0 as hinted at by the brightness gra- 
dient. The emitting gas density can be determined from 
the fitted emission measure and a simple geometric esti- 
mate for the volume. The east rim spectrum comes from a 
region ~10" long by ~3" thick. Assuming the depth in the 
line-of-sight is the same as its length, we find a volume of 
6.4 X 10'"'''' cm3, which results in a hydrogen number density 
of 23 cm~3. The southeast rim spectrum yields a similar 
density estimate of 10 cm~3. We estimate the emission 
volume of the west rim spectrum (extracted from a region 
of size 50" X 15") to be ^1.5 x 10'^* cm3 which leads to a 
density estimate of 1 cm~3. There is more than a factor 
of 10 variation in the ambient medium density from east 
to west. We can estimate "shock" ages from the densities 
just derived and the ionization timescales in Table 2.1. We 



find values of 55 yrs for the east rim and 500 yrs for the 
southeast rim, albeit with large errors. This suggests that 
the SN shock wave has encountered some portions of the 
bright eastern limb more recently than others. To conclude 
this discussion of the remnant's environment we note that 
the open cluster HS114 (Hodge & Sexton 1966) hes on the 
high density side of the remnant; the bright X-ray limb 
(near the east rim region) is only some 8" (2 pc) from the 
edge of the cluster. 

The complex morphology of SNR 0506—68.0 revealed 
by our new Chandra data makes it clear that estimates 
of its dynamical state based on global properties (e.g., 
size, brightness, integrated spectrum) arc questionable. 
We make the unjustified but plausible assumption that 
the remnant's evolution toward the low density northwest 
hemisphere can be described by an abiabatic Sedov-phase 
model. The radius of the remnant in this direction is ~50" 
(12 pc) and the ambient density, using the results in the 
previous paragraph, is ~0.25 cm~3. With these values and 
a canonical explosion energy of 10^"^ ergs the age of SNR 
0506-68.0 is -4600 yrs. 

The high density environment, proximity to a young 
star cluster, and presence of enhanced low-Z metal abun- 
dances in the remnant all point toward a core collapse SN 
and massive star progenitor for SNR 0506—68.0. A neu- 
tron star (or, more speculatively, a black hole) is therefore 
a plausible identification for the compact central X-ray 
star. The lack of an extended nebula and the relatively low 
radio-to-X-ray flux ratio argue against a young, rapidly- 
rotating, high spin-down-rate pulsar. However the X-ray 
luminosity (a few times 1033 erg s~^) and spectrum (e.g., 
kTBB ~ 0.75 — 1.3 keV) are broadly consistent with the 
emerging class of central compact X-ray emitting objects 
such as the one in Cas A (Pavlov et al. 2000) that are 
believed to be isolated neutron stars. Unfortunately the 
ciirrcnt optical magnitude limits cannot definitively rule 
out an extragalactic origin for the X-ray source, although 
the star's location very near the center of SNR 0506—68.0 
gives strong support for a physical association. An un- 
ambiguous identification of the X-ray source will require 
further study across the electromagnetic spectrum. 

This paper utilizes public domain data originally ob- 
tained by the MACHO Project, whose work was per- 
formed under the joint auspices of the U.S. Depart- 
ment of Energy, National Nuclear Security Administra- 
tion by the University of California, Lawrence Livermore 
National Laboratory under contract No. W-7405-Eng-48, 
the National Science Foundation through the Center for 
Particle Astrophysics of the University of California un- 
der cooperative agreement AST-8809616, and the Mount 
Stromlo and Siding Spring Observatory, part of the Aus- 
tralian National University. We also acknowledge J. Dickel 
and P. Ghavamian, who helped with various aspects of 
this research. Partial support was provided by Chandra 
grant no. GO2-3069X to Rutgers University, NASA grant 
no. NAG5-9281 to SAO, and NASA Contract NAS8-39073 
to SAO. 
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